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Методика выделения топологических 
признаков на трехмерных изображениях 


Представлена новая методика получения и анализа топологических характеристик для последующей 
классификации на трехмерных изображениях. Для получения топологических структур использовались 
трехмерные дистанционные карты, анализ применялся с помощью рассмотрения совместных свойств 
скелета, границы и выпуклой оболочки объекта. 


Введение 


Трехмерные объекты на изображениях отличаются большой вариабельностью форм 
и структур, характеристика которых является сложной задачей. Характеристика структур 
является очень важной проблемой в задачах исследования материалов, мониторинга за- 
болеваний, контроля качества. 

Выделение средних линий объектов (скелетизация) является важным этапом в таких 
задачах для последующих вычислений различных (геометрических, топологических) 
характеристик объекта, а также для описания и распознавания объектов. В медицинских 
приложениях скелетизация может использоваться для мониторинга заболеваний путем 
контроля изменений на изображениях. 

В большом количестве работ, посвященных проблеме выделения топологических 
признаков, рассматриваются такие топологические примитивы: окрестность, связность, 
путь, полость и т.д. [1]. Они, безусловно, важны, однако для достаточно сложных 
объектов получается слишком большая дескрипторная часть, что осложняет анализ 
объектов. 

В нашей статье рассматриваются более сложные топологические характеристики, 
которые позволяют достаточно полно описать особенности сложных ветвящихся объектов. 

В основе вычисления использовался скелет, который имеет следующие основные 
свойства: 

— сохраняет топологию оригинального объекта; 

— расположен в центральных точках геометрической области объекта; 

— все воксели скелета связаны друг с другом; 

— сохраняется после вращения изображения (по крайней мере для угла, который 
является кратным д/2); 

— имеет ширину равную 1 вокселю; 

— для каждого вокселя скелета известно расстояние до границы объекта. 

Последнее свойство основано на построении скелета с помощью дистанционных 
карт и позволяет эффективно вычислить топологические характеристики объекта, как 
будет показано далее. 

Целью данной работы является разработка методики получения и анализа то- 
пологических характеристик с помощью рассмотрения совместных свойств скелета, 
границы и выпуклой оболочки объекта. 
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Трехмерные дистанционные карты 


Рассмотрим п-мерное бинарное изображение /[, для которого определены мно- 
жества вокселей (элементов объемного изображения) объекта О и вокселей фона В. 
Пусть хе №" - некоторая координата, тогда через /(х) обозначим значение вокселя 
с координатой х в изображении /. 

Дистанционное преобразование — преобразование, которое каждой точке изображе- 
ния ставит в соответствие расстояние в заданной метрике до ближайшей точки фона [2]. 
Формальная запись: 


0, если Г(х) Е В 


91 у (1) 
‚УЦКу)Е В, если [(х) ЕО 


Та(х) = 


ших —у 


где | — в — некоторая метрика. 


Также существует понятие знакового дистанционного преобразования. В этом 
случае дополнительно вокселям фона ставится в соответствие число, противоположное 
расстоянию до ближайшей точки объекта. То есть: 


=ши(к-у ‚УКу) Е О, если Г(х) Е В 
‚УКу)Е В,если [(х) ЕО 


14(х)= (2) 


ших —У 


Точки хребтов на дистанционной карте соответствуют вокселям, которые 
находятся в центре объекта. Они выступают в качестве потенциальных кандидатов 
точек скелета. Ниже перечислены некоторые подходы, использующиеся для поиска 
вокселей-кандидатов: 

—методы утоньшения, использующие дистанционную карту для определения 
приоритета выбора вокселей для удаления [3]; 

—методы поиска градиента [4], предполагающие выявление районов неодно- 
родного градиента и маркировки таких вокселей в качестве кандидатов на удаление; 

—методы вычисления дивергенции используют в [5] в качестве функции прио- 
ритета для удаления «простых» вокселей с малым значением дивергенции; 

—методы адаптивного утоньшения характеризуются сравнениями между значе- 
нием дистанционной карты в вокселе и средним значением дистанционной карты его 
соседей [6]. 

Множество вокселей-кандидатов обычно имеет большую размерность и следу- 
ющий шаг, как правило, заключается в уменьшении их количества. 

Для связности большинство алгоритмов используют минимальные остовные 
деревья [7], кратчайшие пути [9] или другие алгоритмы на графах. Некоторые методы 
сначала используют объединение, затем удаление вокселей путем нахождения крат- 
чайшего пути в связанном множестве [4]. 

Основное преимущество этих методов в том, что вычисление дистанционной 
карты происходит очень быстро (за линейное время от числа вокселей), что очень 
важно во многих приложениях. 


Выделение топологических характеристик 


Топологические характеристики описывают форму объектов и их морфологи- 
ческую взаимосвязь. Наиболее яркими и интересными из них являются следующие: 
число Эйлера, характеризующее вложенность объектов внутри измеряемого; выпуклый 
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периметр, соответствующий периметру минимального выпуклого полигона, в который 
вписывается объект; количество дыр или отверстий внутри объекта. 

Покажем, как можно эффективно вычислить топологические характеристики с 
помощью построенных скелета и контуров. 

Иногда встречаются объекты, которые чувствительны только к определенному 
типу шаровидности, который свидетельствует о локальной нерегулярности и вычисляется 
как отношение объема объекта к объему эквивалентного шара: 


Зх уомте 


&1оБо5йу = (4) 


4х лх га из ` 
Здесь радиус можно посчитать как половину максимума расстояния между вок- 
селями скелета. Радиус можно вычислить за линейное время, сравнивая все воксели 
сегмента скелета с концевыми вокселями. 
Еще одной характеристикой, описывающей форму, является шарообразность, 
которая определяется как отношение радиусов описанной и вписанной окружностей 


(рис. 1): 


эрйейсйу = = : (5) 
ош 
где А» — внутренний радиус; Кош — внешний. Внутренний радиус берется как 
максимум значения толщины сегмента. Обозначим точку, в которой достигается 
максимум толщины, как Риах. Внешний радиус берется как максимум расстояния от 
Ртах до границы объекта. 


Ком 


Рисунок 1 — Внутренний и внешний радиус объекта 


Несмотря на большое количество характеристик, описывающих форму объекта, для 
характеристики выпуклости (сопуехИу) используется отношение обычного объема к 


выпуклому (рис. 2): 


уоште 


(6) 


сопуехйу = — о. 
сопуехуоите 


Если объект выпуклый, эта величина будет равна 1, при появлении впадин она 
уменьшается. Здесь объем представляет собой сумму площадей объекта на каждом 
слое с учетом масштабирующего коэффициента. 

Похожая величина используется при изучении плотности объекта. Называется 
она «целостность» (зоПАИу) и определяется как отношение обычной и выпуклой 
площади поверхности (рис. 2): 

агеа 


5оПайу = (7) 


сопуех агеа 
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Здесь площадь поверхности вычисляется как сумма периметров контуров объекта 
на каждом слое. Если объект целый и не имеет дыр и впадин, эта величина равна 
единице, при их появлении она уменьшается. 


Ё ол 
> 
: 
Е 
д 0.5 
| 
0.1 03 0.5 07 09 
Целостность 


Рисунок 2 — Объекты разной выпуклости и целостности 


Когда анализируются изображения со многими объектами, например, параметры 
формы вычисляются индивидуально, описание общей формы выполняется на основе 
индивидуальных параметров, например таких, как стандартное отклонение площади, 
основная и дополнительная оси, периметр и округлость или среднее компактности, 
элонгации, периметра и площади. Все эти величины предназначены в основном для 
описания площадных объектов типа клеток. 

Для протяженных объектов, например сосудов и волокон, в продольном сечении 
многие измерения имеют свою специфику (рис. 3), например длина обычного объекта 
определяется как максимальное расстояние между двумя точками объекта, в то же 
время основной физический смысл имеют длина и толщина развернутого, вытянутого 
волокна. Такая длина соответствует размерам средней линии, а толщина — минимальному 
расстоянию между двумя граничными точками внутри объекта и может определяться как 


ИБгейепай = зкееюоп еп ай. (8) 


Длина скелета подсчитывается следующим образом. Перебираются последова- 
тельно все воксели скелета с одного конца до другого. Пусть (4х, 4у, 42) — смещение 
текущего вокселя от предыдущего. 

Следующая таблица описывает, насколько изменяется текущая длина скелета, 
в зависимости от расположения соседних вокселей скелета. 


Таблица 1 — Процедура подсчета изменений текущей длины скелета. 


Условие взаимного расположения Изменение текущей длины скелета 
аз(ах) + аб5(4у) + аб5(42) == +1 
аБ5(Ах) + аБ5(4у) + аБ5(42) == 2 о, 
аз (ах) + аБ5(4у) + аБ5(42) == +5. 


В отличие от двухмерного случая, где толщина волокна является частным от 
деления площади на длину, в трехмерном случае из-за дискретности между слоями 
можно говорить о максимальной, минимальной и средней толщине волокна: 


тах йпит дге ил = тах {у | у Е 5Кееюп} ; (9) 
ши йиит ПБге илй = пота {у | у Е зКееюп} ; (10) 
ауегазе ПЬге ила = ас {у | у е зКе[еюп}. (11) 
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Толщина волокна, 


Длина объекта, 


Длина волокна 


Толщина объекта 


Рисунок 3 — Отличие между линейными характеристиками объекта и волокна 


На основе разницы между этими понятиями определяют ряд новых параметров, 
например закрученность (си!) волокна — отношение длины объекта к длине волокна: 
еп2й 


и =—. (12) 
ИБге [еп й 


Здесь ДБге [еп определяется по формуле (8), а в качестве [епей берется 
расстояние между концами сегмента скелета. 

Помимо перечисленных, еще одной топологической характеристикой является 
элонгация (е]опгайоп), которая определяется как отношение длины и ширины объекта 
(рис. 4): 

УЛ 


, 13 
1епа1й и 


еопгайоп = 
где ий — ширина описанного прямоугольника; /еи21й — его длина. 
Преимуществом использования скелета вместо исходного объекта является 
меньшее количество операций, необходимое для подсчета ширины и длины прямо- 
угольника. 


Величина этого параметра меняется от 0 до 1, если она равна 1, то объект имеет 
квадратную или круглую форму. 


а) 6) в) 


Рисунок 4 Объекты с разной элонгацией: 
а) высокая элонгация; 6) и в) низкая элонгация 
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Данное определение элонгации не может качественно описать форму искривлен- 
ных объектов, для которых длина не соответствует основной оси. В таком случае 
элонгация определяется как отношение объема объекта к кубу его настоящей толщины, 
которая соответствует максимальной области утоньшения: 


6 хуомте 


ЕЕ ь (14) 


еопзайоп = 


где 4 — максимальное значение толщины на рассматриваемом сегменте скелета. 

На практике большой информативностью обладают комбинации параметров 
между различными объектами, являющиеся составляющими одного целого, так называ- 
емые межобъектные характеристики. Для описания таких комбинаций можно определить 
следующие характеристики: ветвистость (Бгапсбтезз), закрученность (си тезз), хвоста- 
тость (а1пе$$), закрученность хвостов (а ситез$$), нормированная хвостатость 
({а15 гано). 

Ветвистость и закрученность описывают сложность формы скелета. Они опре- 
деляются как отношение числа узлов или сегментов скелета к длине скелета соответ- 
ственно. Сегменты скелета определяются как фрагменты, ограниченные узлами или 
точками окончания скелета: 


Бгапсйтез$ = ОИ $ (15) 
зкееоп (еп 2 й 

си [тез5 = ровное : (16) 
соипК5еетепт5) 


Определение хвостов позволяет определить характеристики для описания зави- 
симости особенностей окончания объекта к его сложности. Дело в том, что в реальной 
жизни объект может иметь не только отдельные ветви, но и пересекающиеся структуры, 
образующие замкнутую сеть. Для определения таких особенностей и вводятся опре- 
деления хвостатости, закрученности хвостов и нормализованной хвостатости. Они 
определяются как: 


о (17) 
зкееоп (еп й 
аси" тез = О ба ; (18) 
соип(5еетепт5) 
1айятано = А з (19) 
зкееоп [еп й 


При этом хвост определяется только своими концевыми вокселями. Для того, 
чтобы сегмент был хвостом, необходимо, чтобы один концевой воксель имел не больше 
одного соседа в 26-окрестности. Длина хвоста из формулы (19) подсчитывается по 
тем же правилам, что и длина скелета из формулы (8). 


Выводы 


Предложенная методика позволяет на основе скелета вычислить топологи- 
ческие характеристики структуры объекта. Широкий набор предлагаемых характеристик 
позволяет не только описать сложность объекта, но и в случае динамических исследо- 


212 «Искусственный интеллект» 4’2011 


Методика выделения топологических признаков на трехмерных изображениях АН 


ваний определять характеристики роста и изменения объекта в зависимости от различных 
физических воздействий, например по ветвистости можно судить о равномерности 
плотности ткани или материала, в котором растет структура. Это особенно важно при 
изучении различных коммуникационных особенностей (снабжение кровью опухоли, 
изменение электрической активности нервной ткани, влияние магнитных полей на 
доменные структуры и т.д.). Данная методика позволяет проводить мониторинг роста 
структур на ЗО изображениях. 

Полученные результаты будут способствовать развитию диагностического про- 
граммного обеспечения лучевых методов для правильного расчета данных по динамике 
регрессии опухолевых образований, оптимизации принятия решения по вопросам 
лечения. Определение правильного пути лечения позволит значительно снизить сто- 
имость лекарств и лечебных процедур. Работа направлена на улучшение информаци- 
онного содержания данных первичного диагноза, а также для обеспечения комфорта 
специалистам-медикам в ходе анализа результатов обследования. Работа частично 
сопровождалась проектами МНТЦ В-1682 и В-1636. 
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А.М. Недзьведь, Д.А. Гончаров, С.В. Абламейко 

Методика вид1лення тополойчних ознак на тривимрних зображеннях 

Представлена нова методика отримання та аналзу тополопчних характеристик для подальшо! класификаци 
на тривимрних зображеннях. Для отримання тополопчних структур використовувалися тривимрн! 
дистанщйн! карти, анал1з застосовувався з допомогою розгляду сшльних властивостей скелета, меж! 
та опукло1 оболонки об’екту. 


А.М. М№е@суеа, О.А. Сопсйаготу, 5.7. АМатеуко 

Тесвшадие Гог зерагайоп 01 ф0оро1ор1са| Ёеабиге$ оп Фе @гее-Чптеп$0па| ппасез 

ТБе рарег 15 4еуое4 {10 Фе рго Мет оЁ оМашите ап4 апа]у71те Фе юро|оз1са[ Ёеабиге$ Рог 1е забзеацет 
с1а55ШсаНоп оЁ Шгее-Читепз1опа] ппазез. То оМашт Ше ‘юроозлса| згасвгез, уе изе4 йгее-4ппеп1опа] 
Чз{апсе тарз, Ше апа]уз1$ аррПе4 Бу сопз14егте е ] ош ргорегиез оЁ фе зКееюп, фе Боипдагу апа ше 
сопуех Ви! оРфе оБуеси. 


Статья поступила в редакцию 22.06.2011. 
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